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Abstract. I present a broad overview of modelling of the Narrow Line Re- 
gion (NLR) of active galaxies, and discuss some of the more recent models we 
currently have for the emission from the NLR. I show why the emission line 
ratios from the NLR are constrained to certain observed values, and describe 
what physical parameters we can derive from observations using emission line 
models. Also presented are some examples of this, looking at the metallicity and 
excitation mechanism of active galaxies. As a final point, the limitations of the 
current models are discussed, and how how the combination of modelling and 
theory can help us solve some of the questions that still remain within the NLR. 



1. Introduction 

The Narrow Line Region (NLR) is the region of extended interstellar gas ionized 
and heated by the active galactic nucleus (AGN). This region is classified as 
"narrow-line" as the gas lies outside the dominating influence of the central 
black hole, and the observed line velocity widths are much less than found in 
typical Broad Line Regions (discussed earlier in these proceedings). Even so, 
typical NLR velocity widths are in the range 200 <FWHM< 500 km s _1 , and 
even higher velocity tails are observed. 

Another distinguishing feature of the NLR is the gas density. In the NLR, 
the gas density is low enough that emission lines arising from magnetic dipole 
transitions can occur. These forbidden lines, like [Om]A5007A and [Nn]A6584A, 
are strong gas coolants and dominate the NLR emission line spectrum. It is 
these lines that are typically used to identify both the NLR and AGN in optical 
spectra. 

The NLR is also much more spatially extended than the BLR, reaching 
kiloparsec scales if the "Extended Narrow Line Region" (ENLR) is included. 
At such distances the NLR extends well beyond the nuclear obscuring material 
commonly called the torus, meaning that it is generally still visible even when 
other more nuclear regions, like the BLR, are hidden. 
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Figure 1. BPT diagram of [0 ill] A5007/H/3 vs. [N n] A6584/Ha, showing the 
distribution of SDSS DR4 emission line galaxies. The three classes of emission 
line galaxies, Star-forming, LINERs and Seyferts, are seen as three separate 
branches, with the empirical divisions marked to guide the eye. The LINERs 
and Se yferts are often classified together as an "AGN branch" JG roves et al.l 
I2006U Figure 1 from ). 



2. Line Ratio Diagrams 



In modelling the NLR, we want to reproduce and explain the geometry and 
spatial extent of the region and, most importantly, the emission from this region. 
The emission lines from the NLR are not only indicative of the region, but also 
provide diagnostics for the density, temperature, and ionization mechanism of 
the NLR gas. 

A useful visual tool for comparing models with observations of many NLR 
are Line Ratio Diagrams. As demonstrated by figure HJ these plot ratios of 
emission lines against each other and, depending upon the emission lines chosen, 
can show clear relationships w ith density, metallicity, and ionization mechanism 

(|Veilleux fc Osterbrocklll987l h Figure [1] shows the "BPT " diagram, one of the 

earlies t and strongest emission line diagrams, suggested bv lBaldwin. Phillips. Terlevicii 
(1981). Thisdia gram is able to distinguish three different classes of narrow emis- 
sion line galaxies: star-forming galaxies, Low Ionization Narrow Emission line 
Region or LINER galaxies, and the Seyfert galaxies which show typical AGN 
NLR emission line ratios. These three classes differ in their ionizing source and 
hence nebular temperature and emission. In terms of the NLR models, the final 
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Figure 2. BPT diagram showing four different models as labelled. Four 
nearby AGN are marked; 3 typi cal Seyferts and M87, known to ha ve a strong 
jet and shocks. (Figure 1 from lAllen. Dopita. &: Tsvetanovlll998t ) 



reproduced spectra must lie in the top right section of this diagram in the region 
occupied by the "pure" Seyfert galaxies. 



3. Emission Line Region Modelling 

To determine the resulting emission lines from ionized nebulae, emission line 
region models need to determine the de nsity, temperature and ioni zation state 
of the ionized gas (discussed in detail in lDopita Sz Sutherland! [20031 ) . The tem- 



perature and ionization are determined fundamentally by the ionization source, 
of which there are two main pathways: Photoionization and Shock ionization. 



3.1. Shock Ionization 

In shock ionization modelling of NLR, the gas is excited collisionally through 
shocks caused by interactions with a jet or winds arising from the central AGN 
source. These models generally take as input the gas abundances and pre-shock 
density, shock velocity V s , and a parameter related to the magnetic field strength 
B/n 1 / 2 . 
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Typical line values for shock models are shown in fig ure [2J These mod- 
els, discussed in detail in I Allen. Dopita. &: Tsvetanovl ( 19981 ). cover the range of 



parameter space expected for Seyfert NLRs. As can be seen from the figure, 
the simple shock models are unable to reproduce typical Seyfert NLR line ra- 
tios. For some LINER-like objects however, they might possibly be the exciting 
mechanism. 



3.2. Photoionization 

In photoionization modelling, the NLR is excited by the UV and X-ray photons 
from a central source, thought to be the accretion disk surrounding the central 
black hole in AGN. These mo dels have b e en ex plored for many years, with 
good reviews being the book by lOsterbrock! (jl989l ) and the lecture by H. Netzer 
(|Blandford et alJll99d in ). 



The typical parameters for photoionization models are the gas abundances 
and initial density, the size or column depth of the model cloud, the incident 
ionizing spectrum, and the incident ionizing flux of the radiation. The final 
parameter is often described in terms of the ionization parameter, U, a dimen- 
sionless measure of the density of ionizing photons over the gas density, 

1 f°° F 

U = — ^du. (1) 

U is a controlling parameter of the emission line spectrum in photoionization 
models. 

The resulting models for a range of typical NLR values can be seen in figure 
[2j The models do reproduce the ratios quite well for a limited range of param- 
eters, as well as others such as [O I]A6300/H<x However these models cannot 
reproduce both low-ionization and high-ionization line strengths simultaneously, 
for example failing to reproduce He IIA4686/H /3 with the sam e parameter set, 
and these simple models haven been ruled out ( Stasinskalll984l ). 



With both simple shock models and simple photoionization models not 
surprisingly ruled out, several more complicated models have been put forward 
to try and better reproduce and explain the NLR. 

3.3. Shock+Precursor 

It was realized early on in shock modelling that fast shocks (V x > 1 50 km s" 1 ) 
would produce ionizing photons ( Dopita Sz Sutherland! 11995 . 1 19961 ) . As post- 



shock gas cools it produces ionizing photons which diffuse up- and down-stream, 
and ionize the pre-shock gas. Thus the ionization in fast shocks is actually a 
combination of both shock and photo-ionization, where the photoionization is 
determined by the shock velocity. 

These shock -l-precursor models are able to reproduce the observations quite 
well, as shown in figure [2j and can also produce some of the higher ioniza- 
tion lines. However, the problem with these models is that they require shocks 
throughout the NLR, meaning shock signatures should always be visible. So by 
themselves these models cannot explain all NLR emission. 
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3.4. Multi-Component Photoionization 

The next level of complexity in terms of photoionization models is multi-component 
or multi-cloud photoionization, where the combination of two or more photoion- 
ization models is used to reproduce both the high and low ionization lines in the 
NLR. 

Generally, most models limit themselves to two components to minimize 
the number of free parameters. These models have been used to ex amine spe- 



cific galaxies (e.g. Kraemer Sz Crenshaw 20001 ; Morganti et al. 1991) and more 



generally trying to explain s pecific line ratios (e.g. Komossa & Sc hulzl 119971 ; 



Muravama &: Taniguchi| [l998) or line strengths ( Baskin Laorl l2005). 

The main problem with these models is that as you increase the number of 
clouds, the problem becomes unconstrained. A way around this is to provide 
a physical basis for the multiple clouds, and hence physical constraints. Here I 
mention three current models that try to deal with this problem. 



Ayi/i Models 



3.5. 

The A M /i models of iBinette et al.l ([l996) consider the NLR to be made up of 
hot, highly ionized, matter boundecQ clouds and, lower ionization, ionization 
bounded clouds. In their model, the ionization bounded (I-) clouds see the 
absorbed spectrum from the matter bounded (M-) clouds, and the I-clouds are 
also a higher density. The resulting emission line spectrum from these models 
is the controlled by the ratio of these two components, A M ^. 

Two possible physical pictures of their model are shown in figure [3] (their 
figure 4). The first (figured) is where the clouds are separate, with inner clouds 
being M-clouds, and the outer are I-clouds, while the second (figure [3b) shows 
composite clouds, with an I-cloud core and a M-cloud photoevaporating halo. 

These models reproduce the strong lines quite well as shown in figure [2] by 
the A M /i marked curve, and they also reprod uce very well strong high ionization 
ratios such as He IIA4686/H/3 (see e.g. fig. 7. IBinette et~allll996h 



3.6. Local Optimally Emitting Clouds 

The Locally Optimally-emitting Cloud model of Ferguson et al. ( 19971 ) is an 



extension of a model for the BLR. It uses the fact that the spectra we observe in 
AGN NLRs are likely to be dominated by selection effects: each line we observed 
arises from the clouds in which it is most strongly emitted. As each line emits 
strongest near its critical density, this recreates the linewidth-critical density 
relation observed in some NLR. 



For their model. iFerguson et al.l run a grid of simple photoionization models 
covering a range of densities and radii from the nucleus (connected to incident 
flux). The total line flux is then the integral of these over the NLR cloud 
distribution function, ifi, 



F\me oc 



model 



V, n)ip(r, n)dn dr. 



(2) 



1 Matter bounded clouds have low enough column densities that not all ionizing radiation is 
absorbed. This compares to Ionization bounded clouds where all ionizing radiation is absorbed. 
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Figure 3. Two possible physical geometries for the A M /i model. (Figure 4 
from iBinette et al"1ll996D 



While ip(r,n) can take any form, for simplicity they assume separate power-law 
distributions, ip{r, n) oc r a n^. With this distribution and simple models they ob- 
tain models that reproduce the observations very well, as shown by their figure. 
However again, these models are actually relatively physically unconstrained, 
and importantly they ignore the fact that the NLR is clumpy in nature and 
dusty. 

3.7. Radiation Pressure Dominated Dusty Clouds 

The dusty, radiation pressure dominated cloud models by I Groves et al.l ( 2004al lbl) 



take in the observations that indicate that the NLR clouds are likely dusty and 
clumpy in nature to come up with a physical model that explains why AGN 
NLRs cluster where they do on line ratio diagrams. It uses the fact that at 
high ionization parameter (U) dust dominates the opacity in dusty gas, and 
that in isobaric systems the gas pressure gradient must match the radation 
pressure gradient to realise that radiation press ure on dust will do minate the 
NLR structure at high U. As shown in detail in Dopita et al. ( 20021 ) this leads 



to a self regulatory mechanism for the loca0 ionization parameter and hence for 
the emission lines. These models reproduce the observati ons very well for bot h 
high and low ionization species, as shown in the figures in iGroves et all (l2004bf h 



and converge in the region in line ratio diagrams where the NLR is observed to 
be. 



2 C^iocai is defined using the local absorbed ration field and density. 
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Figure 4. Possible p hysical picture for a NLR cloud in the Dusty, P ra d 
model. (Figure 1 from iDopita et al.ll2003 ) 



The physical picture for these models is shown in figure [U revealing the 
main emitting region, the observed NLR tails, and possible sources for the very 
high ionization "coronal lines" observed in some AGN. In many respects it is 
very similar to the Ay^n picture shown in figure [3)3. 

These models have several other benefits arising from the dust, such as a 
hardening of the radiation field (cf. the A M n models), increased temperature 
through photoelectric heating, and of course the corresponding IR emission. 
This p rovides another constraint of NLR models, as the models (e.g. iGroves et al 



2006a) must match r e cent o bservations of IR emission from the NLR (see e.g. lBock et al. 



(200Q); M ason et all ([2006) and several other papers in this proceedings such as 
by Schweitzer). 

The dust is also a hindrance in some respects, as Iron lines like [Fe VII] tend 
to be too weak due to dust depletion. Similarly, other high ionization coronal 
lines are unable to be reproduced. It should be noted however, that these models 
are for individual NLR clouds, not the NLR as a whole. 



4. Low Metallicity AGN 

One of the uses of NLR models is to determine the average physical quantities 
of a specific AGN. Inversely, it is possible to use NLR models to predict what 
type of emission AGN of a specific type would have. For example, in a recent 
work [Groves et ail ( 2006bl ) used NLR models to determine strong metallicity- 
sensitive emission line ratio diagrams, and where low metallicity AGN would 
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Figure 5. The effects of low metallicity on AGN. Fig. a) shows NLR models 
of decreasing metallicity (as labelled) on the BPT diagram. The underlying 
image shows the distribution of the SDSS DR4 emission line galaxy sample. 
Fig. b) shows how the mass-metallicity correlation applies in AGN as well, 
showing h ow the galaxy mass v aries with position in the BPT diagram (as 
define d inlGroves et al. (|2006b)) which correlates with metallicity. (Figures 
from iGroves et al.ll2006bf ) . 



appear on these. Figure [5^ reveals one possible diagram, the BPT diagram, 
well known to be metallicity se nsitive due to the secondary nature of nitrogen 
( Storchi-Bergmann et al1ll998l ). These diagrams were then used to look for 
possible low metallicity AGN in the SDSS emission line galaxy sample. Only 
40 likely low metallicity Seyferts were found in a sample of ~23000 Seyferts (or 
~ 8800 "pure" Seyferts). 

In figure [SJd) the variation of the gal axy host mass with p osition in the BPT 
Seyfert branch is shown (details found in Groves et al.ll2006bl ). This reveals that 
the mass-metallicity correlation holds in AGN as well (as expected). 



5. Diagnosing Emission Line Galaxies 

Emission line ratio diagrams were originally suggested to separate Star-forming 
galaxies from Seyfert type-2 (NLR only) galaxies, as both show strong nar - 
row lines, although at different strengths (Baldwin. Philli ps^ Terlevichlll98"ll ). 



With the advent of large spectroscopic surveys like SDSS this work has pro- 
gressed much further, with the BPT diagram showing a clear continuum of 
objects with the edges of the V-Shaped figure corresponding to galaxy evolution 
effects, ionization limits, and observation limits. While the distinction between 
the star- forming galaxies and AGN is clear in the BPT diagram (Figure [T]) apart 
from at low [OIII]/H/3, the distinction between the LINERs and Seyfert galax- 
ies is not so c l ear, a nd these objects are often classed together. Recent work by 
iKewlev et al.l (|2006l ) however, has shown that a clear divide does exist between 
Seyfert and LINER objects as shown in figure This figure demonstrates that 
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Figure 6. The separation of LINERs from Seyferts. Figure a) shows the 
AGN branch of the SDSS galaxies on the [S IIJAA6713, 31 line ratio diagram, 
and the placement of the arcs to measure the separation of the Seyfert and 
LINER branches. Figure b) shows the clear separation of the Seyfert and 
LINER branches in the [S II]AA6713, 31 and [OI]A6300 ratio diagrams. 



LINERs and Seyferts appear as distinct populations in the [S II] and [O I] line 
ratio diagrams. 

With this distinction, the separate galaxy properties of all thre e classes have 



been ex amined and compared. In addition, wi th NLR mode l s like [Groves et~al 



( 2004al fbh and Star- forming galaxy models like lDopita et al" ( 20061 ) we are com- 
ing a far way in understanding the full spread of emission line galaxies in line 
ratio diagrams. 



6. Conclusions 

With current Narrow Line Region models we have a reasonably good physical 
understanding of the NLR and its appearance and emission lines. These models 
can be used to determine properties of the NLR and AGN, but caution is still 
needed, as these models are limited. 

An individual NLR consists of many clouds, with a probable range of den- 
sities, pressures and incident spectrum. These clouds could also vary in their 
abundances, metallicity and dust properties, and certainly vary in shape and 
size. The differences between different active galaxies will be even larger. Thus 
models must account for this, but similarly they must account for the fact that 
NLR spectra are all very similar. The current NLR models are almost there but 
a full physical, geometrical model is not yet available. 
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